Temperature and Field Dependence of the Energy Gap of MgB2/Pb planar junction 
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We have constructed MgB2 /Pb planar junctions for both temperature and field dependence stud- 
ies. Our results show that the small gap is a true bulk property of MgB2 superconductor, not 
due to surface effects. The temperature dependence of the energy gap manifests a nearly BCS-like 
behavior. Analysis of the effect of magnetic field on junctions suggests that the energy gap of MgB2 
depends non-linearly on the magnetic field. Moreover, MgB2 has an upper critical field of 15 T, in 
agreement with some reported Hc2 from transport measurements. 
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I. INTRODUCTION 

Recently Nagamatsu et al.0 has discovered supercon- 
ductivity in the commonly available compound MgB2 
with a Tc of 40 K. Similar to the cuprates, MgB2 has 
a layer structure and hence many of its superconducting 
properties may show anisotropic effect. For instance, the 
anisotropy ratio 7 =j-.£ab/^c has a reported value that 
varies from 1.1 to 9.OB0. There is also evidence that the 
energy gap value is very different along these two dire| 
tions showing anisotropic s-wave or two-gap behavioiu" 
On thCpOjther hand, as clearly dempastrated by isotope 
effcctlijO and neutron scatteringE3Ej, MgB2 is different 
from the cuprates and its Cooper pairs are phonon me- 
diated. 

Although the pairing mechanism in MgB2 is thought 
to be phonon mediated, there are still many experimental 
results that lack appropriate explanation. Many of these 
unanswered problems may lead to unexpected and inter- 
esting physics of superconductivity. For example, there 
is no consensus about the magnitude of the energy gap 
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data on MgB2, as in the case of many other newly dis- 
covered superconductors, are obtained from 
junctions likft,a|CapiiU|ig tunneling microscope 
point contactLfcln3~Ea and tunneling junctionst 

It is critical to determine, whether the small gap value 
reported by many groupstJ is a real bulk property or a 
result of surface degradation. One direct method is to 
measure the temperature dependence of the energy gap. 
Since the structure of a mechanical junction will change 
as the temperature is varied or an external field is ap- 
plied, it is not stable enough to study temperature de- 
pendence of the energy gap. The situation will be worse 
if the sample is not homogeneous and the gap value varies 
with the probe position. The more reliable measurement 
for temperature dependence of the energy gap is from 
sandwich type planar junctions where any variation in 
the tunneling spectra will be a pure result of the sam- 
ple under study not due to any structural changes in the 



junction. To our best knowledge, in this paper we re- 
port the first energy gap temperature and magnetic field 
dependence of MgB2 by planar junctions. 



II. EXPERIMENTAL DETAILS 

MgB2 sample is prepared by reacting Mg turnings 
(99.98%) and boron powder (99.99%, -325 mesh) with the 
stoichiometric composition 1:2 respectively. Magnesium 
and boron are mechanically pressed and sealed in a tantu- 
lum tube (99.9%, 2.4 mm inner diameter). The tantalum 
tube is then sealed inside a quartz tube and placed in- 
side a box furnace at 950 °C. It is then quenched to room 
temperature after two hours. The polycrystalline MgB2 
is then characterized using X-ray diffraction, resistivity 
and dc SQUID magnetometer (Quantum Design MPMS) 
measurements . 

Junctions are constructed by attaching two leads to 
MgB2 sample and molding it inside epoxy resin. It is 
then ground to expose the sample and mechanically pol- 
ished to a smoothness of 0.3 micron. Pb, a superconduc- 
tor with Tc w 7.2 K and Hc{0) « 0.08 T, is evaporated 
on the top as a counter electrode. We used Pb to sharpen 
the peak features and also as a control to monitor the tun- 
neling conditions. In this paper we will limit our analysis 
only to the data when Pb is normal, which is simpler to 
understand. We have also attempted to grow artificial 
barrier by sandwiching a thin oxidized aluminum layer 
between the sample and Pb electrode. This will in gen- 
eral lead to very large junction resistance, even with the 
minimal thickness of aluminum layer. So far, the best 
junctions are still from those with natural barrier. The 
junctions show stability against any temperature changes 
in the full range from 4.2 K to room temperature, but can 
survive only up to a magnetic field (perpendicular to the 
barrier) of approximately 3.2 T. 



III. RESULTS AND DISCUSSIONS 



X-ray diffraction pattern shows no trace of other 
rases in the sample. Both resistivity and SQUID mea- 
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FIG. 1: Magnetization divided by an applied filed of 10 
Oe as a function of temperature for botli Zero Field Cooled 
(ZFC) and Field Cooled (FC) modes. The inset shows the 
magnetization curve M{H) measured at T = 5K. For clarity, 
the X-axis (inset) is limited to small values of H. 



surements have determined the T°"''''* to be 39.5 K (de- 
fined by 2% criteria) with a sharp transition width of 
0.7 K (10%-90% criteria). Fig. (main panel) shows the 
temperature dependence of susceptibiUty for both zero- 
field cooled (ZFC) and field cooled (FC) modes. Tak- 
ing in account the demagnetizing factor of the measured 
cylindrical sample with 7=1 (ratio of length to diam- 
eter), the sample shows a perfect diamagnetic shielding 
M/H = — 1. This result along with a Residual Resis- 
tance Ratio RRR = R{300)/R{Tc) = 8 reflects both the 
sample's good quality and grains coupling. The small FC 
susceptibility signal observed here is a common feature 
for such poly cry gjtallipe MgB2 samples sintered around 
950 °C or highercSEJ. This can be attributed to large 
trapping of flux by cradffi and voids that reflects also the 
good coupling of grainsEJ. It is interesting that the lower 
critical field Hci{5K) = 0.2 T as estimated from the mag- 
netization curve (Fig. |l|, inset). This value is signific anjdy . 
larger than those reported by other research groupscJ'EHl 
and can be attributed to the good quality of the sample. 

The inset of Fig. |^ shows how the conductance curves 
evolve with temperatures below lead Tc- For this S/I/S 
junction we can roughly estimate the energy gap of 
MgB2. It is clear that the spectra are sharpened sig- 
nificantly as the Pb gap Apt, opens up. Since Apf,(0) is 
about 1.2 meV and the peak position of the 4.2K curve 
is at 3.2 meV which can be considered as the half sum of 
gaps, we can estimate AMgB2 to be about 2 meV. The 
corresponding value of 2A/kTc is only about 1.18, much 
smaller than the BCS value of weak coupling supercon- 
ductors. This is consistent witt other small gap results 
from tunneling measurcmcntsE3. As can be seen from 
Fig. 1^ (main panel) , the smaller peak around 9 meV sur- 
vives for temperatures up to 21.16 K. Similar features are 
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FIG. 2: Temperature dependence of the experimental tun- 
neling conductance spectra normalized by the conductance 
at 15 mV (dots). The spectra at temperatures above Tc of 
lead are represented in the main panel with the corresponding 
two-gap BTK fittings (red lines). For clarity the curves are 
vertically shifted. 



commonly found in other tunneling dataO. 

It is intriguing to find this very small energy gap in a 
high Tc superconductor like MgB2. All tunneling data 
published so far can be summarized into three npain-cat- 
egories, |according to the interpret|ation: one-gapEj~EElE2, 
two-gapElEl~Ll and gap anisotropyO. If there is only one 
single gap, then it is likely that the small value is a result 
of surface degradation. However, our data do not support 
this explanation because the gap exists up to the bulk Tc- 
So far, there is no direct observation of two distinguished 
gaps in the same tunneling spectrum. Mostly a small 
feature at a higher energy is interpreted as the second 
gapQ. This two-gap approach is in accordance wilh the 
2D and 3D Fermi surfaces proposed by Liu et al.EJ. On 
the other hand, Chen et al.l3 proposed an anisotropic s- 



wave pairing model with A^y = 5 meV and — 8 meV 
to best fit their tunneling curves. |__. 

We have used Blonder, Tinkham, and KlapwijkEil 
(BTK) model to analyze the curves when Pb is normal 
(Fig. 0, main panel). In this model, in addition to quasi- 
particle tunneling, the possibility of Andreev reflection 
and normal reflection by the barrier is also included. As 
indicated by our fitting, the barrier strength Z of our 
junction is not strong enough to prevent Andreev reflec- 
tion from happening. A depairing term F is also included 
because of shortening in quasiparticle lifetime by differ- 
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FIG. 3: Temperature dependence of the two gaps Ai and 
A2 normalized by their values at T = 7.78 K. The absolute 
values are shown in the inset. The solid line represents the 
expected BCS A(T)/A(0) with = 39.5 K. 



ent scattering processes. 

Here we have assumed the two-gap model to cover also 
the small feature around 9 meV (Fig. main panel). 
We consider this feature as the second energy gap A2 in 
MgB2. Also, the two gaps contribute to tunneling in- 
dependently. Therefore, we assign a small percentage of 
tunneling C2 for A2 and Ci = 1 — C2 for the smaller gap 
Ai. The parameters F, Z, Ai, A2, Ci and C2 are used to 
best fit the curve at 7.78 K and their values are 0.95 meV, 
1.33, 1.75 meV, 8.2 meV, 0.94 and 0.06, respectively. All 
these parameters except the A's are kept constant for all 
higher temperature curves, i.e., the A's are the only ad- 
justing parameters. From the quality of the fittings, it 
is justifiable to say that F, Z and C's are independent of 
energy and temperature within the ranges of our mea- 
surement. Furthermore, the zero bias offset is purely a 
result of Andreev reflection at the barrier. 

The inset of Fig. ^ shows the temperature dependence 
of the two gaps. The main panel shows both gaps nor- 
malized to their values at T = 7.78 K along with the BCS 
expected behavior (solid line). As can be seen, both Ai 
and A2 survive to Tc of bulk MgB2 with a small devi- 
ation from BCS as we try to apply the two-gap model. 
Since the tunneling features are mainly due to Ai (Ci = 
0.94) and they survive up to Tc of MgB2, it is justifi- 
able to consider Ai as a true bulk property of this su- 
perconductor. This two-gap fitting gives A2(Q)iAi(0) w 
4.5, close to both the iteoretically predictedEj and ex- 
perimentally suggestedlij value. Nevertheless, there are 
still unexplained problems with this two-gap model. For 
example, why the large gap contributes that little to tun- 
neling, C2 = 0.06, for such a polycrystalline sample? Our 
analysis above still holds for a single gap superconductor 
by setting C2 = 0, but then we have to explain why Ai 
is so small. 

To further characterize the junctions, we have also 
studied the field dependence of the tunneling spectra 
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FIG. 4: Magnetic field dependence of the experimental tun- 
neling conductance spectra normalized by the conductance at 
15 V. The inset shows the field dependence of Ai measured 
at 4.2 K. 

at 4.2 K (Fig. ^ main panel). The junction in an ex- 
ternal field normal to the barrier is not as stable as 
its performance against temperature changes. It expe- 
riences slight changes even when a small field is ap- 
plied. This can be seen from the development of the 
zero bias conductance peak, similar to that observed by 
another grouped. This can be explained by enhance- 
ment of micro-shorts through the barrier as a result of 
the applied field. Furthermore, most junctions collapse 
and the tunneling spectra transit from quasiparticle to 
Josephson tunneling at fields of about 3.2 T. In this pa- 
per, we focus only on quasiparticle tunneling spectra. It 
can be seen from Fig. |j (main panel) that the quality of 
the spectra has severely degraded when H exceeded He 
of lead. The curve at 0.43 T is more severely smeared 
and depressed as compared to the curve at T = 7.78 K 
(Fig. H). This reflects the fact that this field is already 
greater than He of Pb. Using the peak position of the 
small gap Ai, we can roughly estimate its dependence 
on H as shown on Fig. H (inset). It is worth noting that 
Ai(OT) - Ai(0.43T) > Apb. This supports the above ar- 
gument on the condition of the curve at 0.43 T. Moreover, 
the energy gap depends non-linearly on the the magnetic 
field. Further work should be done to investigate this 
dependence. 

Since MgB2 is a type II superconductor, the effect of 
magnetic field is to produce vortices and hence the order 
parameter on the surface is not homogeneous anymore. 
In a simple model, the tunneling spectrum is an ensem- 
ble of all different gap values sampled within the junc- 
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FIG. 5: Magnetic field dependence of the minimum conduc- 
tance. The finear fit intersects the normal conductance line 
in a point corresponding to Hc2 ~ 15 T. The intersection with 
the vertical axis matches the minimum conductance offset of 
the spectrum at 7.78 K and T (inset). 



tion area (0.05 mm^). If we consider the vortex core as 
a normal region, and the number of vortices produced is 
proportional to the applied normal field, then the zecp 
bias conductance should be proportional to that fieldE3. 
To study this dependence, we have fitted our tunneling 
spectra within the gap by parabola to remove the zero 
conductance peak. The zero conductance can then be 
estimated from the parabola. In Fig. || (main panel) we 
have plotted the zero bias conductance versus the exter- 



nal applied field. It is clear that the zero bias offset in- 
creases linearly with the external field. By extrapolation 
to the normal conductance, we can estimate Hc2 of MgB2 
to be about 15 T. This value is in agreement with Hc2 of 
bulk MgB2 from transport measurements (see, e.g., Ref. 
29, 36) rather than the small reported value of about 
6 T from tunneling analysis (see, e.g., Ref. 24). From 
Fig. H we can also estimate the S/I/N zero bias offset at 
zero field to be around 2.1 mS. This agrees with the zero 
bias offset in the zero field conductance curve at 7.78 K 
(Fig. |, inset). 



IV. CONCLUSION 

We have prepared MgB2/Pb planar junctions to study 
the temperature and field dependence of the energy gap 
of MgB2. The temperature dependence data indicate 
that the small energy gap we have measured is indeed a 
bulk property of MgB2. Moreover, our data do not con- 
tradict the two-gap scenario by considering our reported 
gap of about 2.0 meV as the small gap. Analysis of the 
effect of magnetic field on the junctions shows that MgB2 
has an upper critical field of about 15 T which is consis- 
tent with most transport measurements of Hc2. More- 
over, the energy gap shows a non-linear dependence on 
the magnetic field applied perpendicular to the barrier. 

This work is supported by NSF Grant No. 
DMR997201. 
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